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Abstract 

A simplified theory of the telescopic oscillations in multiwalled carbon nanotubes 
is developed. The explicit expressions for the telescopic force constants (longitu- 
dinal rigidity) and frequencies for telescopic oscillations are derived depending 
on interacting tubes diameters and lengths. The contribution of small- amplitude 
telescopic oscillations to low temperature characteristics of nanotubes such as 
specific heat are estimated. 

1. Introduction 

Multi- walled carbon nanotubes (MWCNTs) are the first discovered nanoscopic 
quasi-lD nanostructures (1). Each MWCNT consists of some nested single- 
walled nanotubes (shells) held mostly by van der Waals forces ((2). 

The telescopic motion ability of inner shells (3) and their unique mechanical 
properties (Q) permit to use multi-walled nanotubes as main movable arms in 
coming nanomechanical devices. The variety of gadgets of this kind was already 
suggested such as a possible mechanical gigahertz oscillator (linear bearing) 
(si sl)> nanoswitch (6!), nanorelay, nanogear j?]), nanorail, reciprocating nano- 



engme 



Therefore the analysis of mechanic characteristics of MWCNT is an 



important objective of study. The present work is just devoted to a simplified 
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continuum version of this problem. The continuum model for telescopic oscilla- 
tions, in which each shell of MWCNT is considered as continuous infinitcsimally 
thin cylinder is described in the next section. The third section devoted to the 
description of the small (thermal) and large-amplitude oscillations for DWCNT 
and MWCNT in the framework of proposed model. Note that the similar con- 
tinuum model was used recently for the investigation of the suction energy and 



large amplitude telescopic oscillations in DWCNT (|9j; llOT ). The contribution 
of temperature-induced oscillations into the tubes heat capacity within Debye 
model is also discussed. In the last section the obtained results are compared 
with the available experimental data (|3j). 

2. Intertube interaction in MWCNT within continuum model 

The interaction energy of two shells of the multi-walled tube is modelled as 
the sum of pair interaction potentials of atoms from different shells. In doing 
so we took for the potential energy of two atoms at the distance I the Lennard- 
Jones potential 

/? m - 76 jl 712 

with attractive and repulsive constants je = 2.43 x 10~ 24 J-nm 6 and 712 = 
3.859 x 10~ 27 J-nm 12 borrowed from (|2j). In accordance with this approximation 
the total intertube interaction energy has form 
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where r^j and T2,j are radius-vectors of the inner and outer tube's atoms re- 
spectively. 

As in (Q) we used instead of (fTJ) the continuum model, for which 
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where r\ and r 2 are inner and outer tubes radii, L\ and L 2 are their lengthes 
(from now on we assume that Li < L 2 ) and Az is the distance between tubes 
outer edges and a is the surface density of carbon atoms in graphene, which 
almost independent on the tube chirality, 

4 _, 
a = — 7= — = 38.2 nm , 
3V3 b 2 

where b = 0.142 nm is the interatomic distance in graphene. 

The integration over variables Z\ and z 2 can be easily carried out analytically, 
but obtained expressions are too cumbersome to be presented here. 

It's trivial, that system energy is minimal when the inner tube is completely 
retracted into the outer tube. In terms of hypergeometric functions the mini- 
mum interaction energy is given by expression 
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Figure 1: The intertube interaction energy for the (5,5)®(10,10) DWCNT with lOnm and 
20nm lengthes 
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3. Telescopic oscillations in DWCNT 



If the outer tube is rigidly mounted, then the longitudinal motion of the 
internal tube is described by primitive equation: 

1 dU(Az) 



a z {Az) 
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m dz ' 

where a is the acceleration of the inner tube with mass m. 

We ignore here the contribution of some defect-induced dissipative forces 
since for high-quality nanotubes they are several orders lower than the retraction 
force due to self-healing mechanism (Q; [ill ) . 

By Q the motion of inner tube is cyclic with period 
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where Azq is determined by relations U(Azq) = U{Li — L\ — Azq) = Eq. 

Due to the special form of potential ([2]) we can separate out two limit form 
of motion: 

1) steady movement for Azo S> b while the potential is just linear in Az; 

2) small oscillations when Azq — > and the potential is just quadratic in 

Az. 

3.1. Large- amplitude oscillations in DWCNT 

In the case of displacement Azo greater than few nanometers the poten- 
tial energy is mostly linear with respect to displacement Az except small- 
displacement zone with quadratic potential energy contribution of which can 
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Figure 2: Longitudinal intertube interaction force for (5,5)@(10,10) DWCNT with lOnm and 
20nm lengthes 
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be neglected. In such a case the oscillatory period can be derived from simple 
formulas for the steady and uniformly accelerated motion. For equal-lengthes 
tubes the period has form 



. 2{E -U min ) v/2to(£o - U„ 
t(E ) = 4W — - = 4- 



where F z is the longitudinal component of retraction force: 

„„ / p | 1 11 i 4rir 2 \ jp I 1 5 i 4r±r 2 

37T 3 a 2 / 21 2*1 I 2'T' 1 ' (r 1+ r 2 )^ J 2*1 I 5 . 2' 1 ' (n+r 2 )a . 

^ = ^ rir2 l v 32 712 76 I " — 

where F z0 is approximately constant for tubes of large radius. 
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Figure 3: Dependence of F z o on tube's radius for the 50nm-lcngth DWCNT 

In terms of the maximal displacement Azq the period can be written as 
follows 

If L 2 > L\ then the steady motion contribution is included and 
,„s y/2m(E - U min ) ( F z L^-Ly 
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Since the period is proportional to the square root of inner tube mass, then 
for the same maximal displacements and the same immobile outer shell the 
oscillatory period of inner tube is proportional to the inverse square of its mass. 
At the same time the oscillatory period does not depend on tubes diameters. 

If the outer tube is also mobile, then the above expressions for periods remain 
valid with mi replaced by the reduced mass 



Note that the interaction energy of atoms forming the tubes rapidly decreases 
with the interatomic distance. Therefore it is enough to consider only interaction 
of adjacent tubes in MWCNT. Since some adjacent shells of MWCNT can be 
rigidly glued by defects, then glued tubes should be considered as double-sided 
shells with integrated masses. 

3.2. Thermal oscillations of DWCNT 

For low temperatures the telescopic oscillations are the least frequency ID 
modes in DWCNT. Therefore for T — > by Boltzmann theorem their energy is 



For small oscillations the maximal potential energy of DWCNT coincides 
with E: 



where k is the rigidity parameter. Taking into account that in the harmonic 
approximation rigidity is the potential energy second derivative with respect 
to the inner tube longitudinal displacement and assuming the tube's radius is 
much smaller of its length for L\ = L 2 yield 
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The harmonic oscillations frequency for this fc(r 1; r 2 ) is 
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and the amplitude of longitudinal thermal oscillations can be estimated as fol- 
lows 



Az = ^M^V6.1(TVr, 



ran. 



We see in particular that Azo is few times lesser than graphene lattice pa- 
rameter. 

To model the intertube interaction force F(n, r 2 , L\, L 2 , Az) in the case of 
L\ ^ L 2 depending on the inner tube edge position Az let us assume that the 
axis of outer tube coincides with the interval [0, L 2 ] of real axis and introduce 
two parameters: 
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where fco is almost independent of n , r 2 for rather large values of these param- 
eters. For n > lOnm we have fcn « —3.7^^-. 

For small maximum retractions of inner tube (Az< 0.3nm) F(n, r 2 , Li, L 2 , Az) 
can be written as follows 



F(r 1 ,r 2 ,L 1 ,L 2 ,Az) = < 



+ kAz, 

0, 



0.3w7i < Az < x 

x < Az < L 2 - L\ - x 



-F Q + k-(L 2 -Li-Az), L 2 -Li-x < Az < L 2 -Li + O.Snm 

(6) 



Here xq = | | is the displacement, which makes the longitudinal retraction 
force equals to zero, k(r\, r 2 , L\, L 2 ) is the DWCNT longitudinal rigidity and 
F (ri,r 2 ,Li,L 2 ) = F(n,r 2 , Li,L 2 ,0). 
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In case of tubes with significantly different lengths (\L\ — L%\ > lnm) the 
formula ([5]) takes on form 

jp ( 1 12 i 4rir 2 | p / 1 6 i 4rir 2 
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and in the case of equal lengths 
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The oscillation cycle can be considered as 
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where Vmaz = y m(-^(^ z o) ~ Umin) is the maximum inner tube velocity and 
U m in = U ( ~l-k2--Li| j j g S y S ^ em minimum potential energy defined by ([3]) . 



With an accuracy of several percent previous equation can by approximated by 



fc \ |Az + a;o| 

Since the system's longitudinal rigidity depends only on the tubes radii and 
inner tube mass is directly proportional to the product of its length and radius, 
then 

IT 

^harmonic ^ y ~J~ ■ 

4. Thermal Oscillation Frequencies in MWCNT 

Considering the long-amplitude oscillations of multi-walled nanotube we as- 
sumed that some nanotube's shells can be bounded by the defects, but in the 
case of thermal oscillations we should take into account the motion of all shells 
because amplitudes of their oscillations are of the same order and much lesser 
then interatomic distance. 

The MWCNT's thermal oscillations frequencies were obtained by solving the 
system of equations for longitudinal displacements of tubes with forces defined 
by the expression above for the intertube potential. In the continuum model 
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arbitrary multi-walled nanotube can be characterized by the inner shell radius 
ro, number of shells n under consideration and the constant distance between 
adjacent shells (d = 0.34n?7j). 

Considering only adjacent shells interaction we find explicit values of the con- 
sequent MWCNT eigenfrequencies LUi,i = l..n. The least eigenvalue is always 
equals to zero corresponding to the whole nanotube translational motion. The 
numerical analysis shows that the maximal frequency depends on the MWCNT 
characteristics but in case of tube with large number of shells (n > 10) it tends 
to the asymptotic value which depends only on the tubes length: 

_ 280 

Umax — UI1Z. 

The obtained value of maximal frequency is seems too low for real tubes because 
defects may rise the longitudinal rigidity of MWCNT. The minimal oscillation 
frequency strongly depends on the number of shells and can be found using the 
following interpolation formula 

4.46 x 10 11 e-° 0128 " ^ 
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Using obtained frequencies we calculated the contribution of tube's telescopic 
oscillations to the individual MWCNT internal energy 



E(T) = V hu t I - + ^ 



where fc^ is the Boltzmann constant and T is an absolute temperature, and 
specific heat 



By © if k B T < hu mm , then 



Hu>n 



while C V {T) ~ nks if ksT S> huj max . 

It is well known that under the Debye temperature the bulk solid specific heat 
decreases as a cubic function of the temperature and for the one dimensional 
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Figure 4: Frequencies of small telescopic oscillations for the 50nm-length MWCNT with 20 
shells 
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Figure 5: The maximal and minimal oscillation frequency for 30nm-length MWCNT for dif- 
ferent number of shells 
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structures such decreasing is given by the linear function. In the case of MWC- 
NTs the telescopic oscillation induced specific heat decreases exponentially in 
the range of small temperatures and the corresponding Debye-like temperature 
is 

where uj m in is given by ([7]). As a result of MWCNTs specific heat exponential 
decreasing their specific heat may be several orders higher that of environment 
for T 3> Tq while for T <C Tb these values change over. For the natural carbon 
nanotubes To varies in the interval 10~ 3 -j- IK. 

5. Summary 

The explicit expressions for longitudinal rigidities and frequencies of small 
and large-amplitude telescopic oscillations of DWCNT and MWCNT were de- 



fl). 



duced in the framework of continuum Lennard- Jones model borrowed from 
The difference between frequencies calculated by obtained formulas and those 
found by numerical methods with account of discrete structure of nanotubes lies 
within the 5%-range. Besides the obtained frequencies of telescopic oscillations 
of MWCNT are in good agreement with available experimental data |a). There- 
fore the considered Lennard- Jones continuum model is seemingly well facilitated 
for description of telescopic trembling of MWCNT. 
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